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Abstract 

Generation of the cosmological baryon asymmetry in SUSY based 
model with broken R-parity and low scale gravity is considered. The 
model allows for a long-life time or even stable proton and observable 
neutron-antineutron oscillations . 



1 Introduction 

Baryogenesis is well known to be a challenging issue in nowadays cosmol- 
ogy and particle physics. We know that, at least within a confidence radius 
of 10 Mpc, or even much more, see pp, the universe is strongly dominated 
by baryons, while practically no antibaryons are observed in the sky. The 
cosmological number density of baryons with respect to the number den- 
sity of photons in cosmic microwave background radiation (CMBR) is quite 
small. From Big Bang Nucleosynthesis (BBN) and the angular fluctuations of 
CMBR, recently accurately measured by WMAP, we know that (see e.g. j2]): 

/3=^^ = ^6-l(T 10 , (1) 

On the other hand, this number is not so small, as it seems, compared to a 
locally baryo-symmetric universe, where the ratio n&/n 7 = ni/n-y would have 
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been some 10 orders of magnitude smaller than this. Explaining the origin of 
the cosmological asymmetry between particles and antiparticles and the value 
of ratio has proven to be an interesting meeting point for particle physics 
and cosmology, especially since, according to the seminal paper by Sakharov 
it requires and gives hints about physics beyond the standard model. 
Indeed, in order for most baryogenesis models to work, baryon number non- 
conservation, C and CP violation, and departure from thermal equilibrium 
are required, thus pushing forward the demand for new fundamental physics. 
In fact even within the context of the minimal standard model (MSM) these 
requirements are fulfilled, but the magnitude of the baryon asymmetry ob- 
tained in the frameworks of MSM is much smaller than the observed value 

There are many extensions of MSM which can explain the origin of the 
cosmological baryon asymmetry and its value, for a review see refs. jU E]. 
Doing that, one has to keep in mind a very strong lower bound on the pro- 
ton life-time j2j: r p > 2 • 10 29 years (model independent) or r p > 10 31 — 10 33 
years depending upon the decay channels. The goal of this work is to present 
a SUSY based model with broken 7?.-parity which allows for a long proton 
life-time, or even for stable proton and can explain the observed cosmologi- 
cal baryon asymmetry. As we see in what follows, it can be done in models 
with low scale gravity where the short distance Planck mass M* is much 
smaller than the standard long distance Planck mass, Mpi = 1.2 • 10 19 GeV, 
due to propagation of gravity in higher dimensions 1 . There is also another 
possibility to have a low Planck mass in the early universe and the normal 
large one at later cosmological epochs, prior to primordial nucleosynthesis, 
due to time variation M* = M*(t). The latter may be induced by the cou- 
pling Rf (0), where R is the curvature scalar and 4> is some time dependent 
scalar field [H|. As a by-product the suggested model may lead to potentially 
observable neutron-antineutron oscillations and B-nonconserving effects in 
processes with heavy quarks. 

If indeed, the effective Planck mass in the early universe is about TeV 
values, we would be in trouble with many versions of baryogenesis scenarios 
which demand much higher energy for their realization, as e.g. GUT or 
baryo-through-lepto genesis. In this case SUSY theories with broken 1Z- 
parity present a viable possibility for successful cosmological baryogenesis 
and, what's more, they allow for some observable effects in particle physics 
with clear signatures and possibly non-negligible magnitudes. 

In the next section the basic picture of the model is given, while the 

^or a string theory inspired, altough non-supersymmetric, realisation of this model 
see ref. 0. 
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third section is devoted to 7£-parity violation. Afterward, in the fourth sec- 
tion, the decay channels of heavy particles which could potentially lead to 
cosmological baryon asymmetry are considered. In section five the out-of- 
equilibrium decays of massive particles in cosmological setup are discussed 
and the magnitude of baryon asymmetry is estimated. Furthermore, the 
processes in particle physics with the variation of the baryonic charge by 
two units, in particular neutron-antineutron oscillations and processes with 
heavy quarks are described (section 6). Finally, the conclusion is devoted to 
a brief summary of the results, and some ideas for future investigations. 



2 Basic features 



In this paper we investigate the possibility of generation of the observed 
baryon asymmetry at relatively low temperatures, around or below the Elec- 
troweak (EW) phase transition, in the context of the Minimal Supersymmet- 
ric Standard Model (MSSM), extended by an addition of 7?.-parity violating 
term in the superpotential, which simultaneously violates baryon-number 
conservation. In this section we give a qualitative picture of how such sce- 
nario could be realised, while in the subsequent sections we go through tech- 
nicalities and details, in particular, we consider a concrete SUSY model with 
B- nonconservat ion . 

Normally the models of low T baryogenesis demand first order phase 
transition in the primeval plasma to break thermal equilibrium because the 
cosmological expansion rate, H = a/ a, at low T is much smaller than the 
reactions rates. Indeed, 



90 / M P , V100/ Mp,' v ; 



where g* is the number of relativistic species contributing into the cosmo- 
logical energy density; at T ~ 100 GeV it is about 100 in MSM and may 
be a factor 2 larger in MSSM. On the other hand, the characteristic rates of 
reactions between elementary particles are: 

T d ~ aNm, and T r ~ 4vra 2 iV 2 T, (3) 

where is the decay rate of a gauge boson with mass m, T r is the typical 
reaction rate, a ~ 0.01 is the gauge coupling constant and N is the number 
of open channels, it could be quite large, N = 10 — 100. 

Evidently, with Mp\ ~ 10 19 GeV, deviations from equilibrium at elec- 
troweak scale for particles with m ~ 100 GeV would be about 10 -17 which is 
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surely not enough to create the baryonic excess (JTJ. We see that in order for 
low temperature models to be efficient, it is necessary to require a lower Plank 
mass, which in turn may be provided by the same theories which present the 
necessary symmetries for keeping the proton stable enough, that is, extra 
dimensions and (super)string theories. Different models of baryogenesis with 
low scale gravity have been considered earlier in ref. j^]. 

Using eqs. (J2J) and (J3J) we find that thermal equilibrium would be strongly 
broken, i.e. H/Y > 1, if the effective Planck mass is bounded from above by 
M* < 10 5 GeV. This estimate is obtained with N ~ f and T ~ m ~ f 00 
GeV. If we demand a weaker condition, H/Y > I0~ 5 , which may be still 
sufficient for generation of f3 ~ I0 -9 , we would see that baryogenesis could 
be possible with M* < I0 10 GeV. 

There could be a complication from another side: if the Hubble param- 
eter is too large, then thermal equilibrium might never be created prior to 
baryogenesis. This may be even favorable for baryogenesis, as it was argued 
in the case of the normal Mpi in ref. fU] , where baryogenesis at reheating 
was considered. A slow process of preheating by the inflaton in the case 
of low effective Planck mass may be supplemented by gravitational parti- 
cle production since the electroweak scale and the particle masses can now 
be close to M*. Correspondingly heavy particles with m ~ M* could be 
efficiently created. Thus we can expect the primeval plasma enriched with 
heavy particles produced by gravity. An analysis of the gravitational particle 
production at the end of inflation can be found in refs. 

The initial properties of the cosmological plasma depend upon the mag- 
nitude of the Higgs field, 0, during the process of (pre)heating. If had 
already reached its nonzero vacuum expectation value, then the produced 
particles would be massive. It is possible that the plasma temperature after 
thermalisation would be below the EW temperature, i.e. the plasma would 
be already in EW broken phase and the EW phase transition never occurred 
in cosmological history. 

In what follows we neglect possible deviations off the standard thermal 
evolution and consider baryogenesis from initially equilibrium state. The 
nonequilibrium effects due to a possible overabundance of heavy particles, 
e.g. t-quarks, would lead to a larger baryon asymmetry. 

Quite efficient breaking of thermal equilibrium in the scenario under con- 
sideration can be created by the following mechanism. At the temperatures 
above the EW phase transition all the Standard Model (SM) particles are 
massless (except for thermal masses, which are in any case lower than the 
temperature at this stage), as well as (possibly) light SUSY particles. At 
a certain stage the phase transition takes place, and the Higgs mechanism 
gives masses to SM particles via spontaneous symmetry breaking. Then, 
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depending upon the final temperature after the transition, the now massive 
particles could get masses which would be higher than the temperature after 
the phase transition, while some other would remain effectively massless with 
m < T. The crucial temperature is known by the standard calculations to 
be related to the Higgs mass, (see e.g. book [12] )• Given the present day 
lower limit on the Higgs mass, we see that the only SM particle which may 
have a mass above the final temperature after the EW phase transition is the 
top quark. We believe that SUSY partners acquire their masses through soft 
symmetry breaking by the analogous Higgs effect at the same EW energy 
scale. 

If the phase transition is fast enough, so that the t-quark would not 
decay down to the Boltzmann suppressed equilibrium density at m t > T in 
the course of the phase transition, then the system would arrive to the state 
when the number density of t-quarks at the end of the EW phase transition 
would be equal to the number density of massless particles and grossly exceed 
their equilibrium value. This makes favorable conditions for baryogenesis 
through t-quark decays which ultimately would reduce their amount down 
to the equilibrium value. 

The rate of the phase transition is determined by the evolution of the 
temperature dependent Higgs mass: 

ml(T) =9^/2 -ml/2 (4) 

where mo is the vacuum Higgs mass and g t ~ 1 is the Yukawa coupling of 
Higgs boson to t-quark. In the course of the phase transition the temperature 
drops down roughly speaking by AT ~ mQ. According to the standard law 
of cosmological temperature red-shift, T = —HT, we find that the time of 
the p.t. is about St ~ H~ 1 {/ST/T) ~ H' 1 . In realist ic case the law of 
temperature variation is different, due to presence of the time varying Higgs 
condensate, but not much. On the other hand, the decay width of t-quark 
is about Y t ~ am t , (see eq. (JIHJ)). Thus the decay rate would be slow in 
comparison with the rate of the p.t. if 

M* < 10 3 m t « 10 5 GeV (5) 

For a bigger M* the t-quark decay would be faster than the rate of the p.t. 
and its number density would follow the equilibrium law, nt ~ exp(— m t /T). 
The baryogenesis through t-quark decays would not be efficient if T/m t <C 1, 
but if this ratio is only mildly smaller than unity then S-nonconserving de- 
cays of t-quark might be efficient enough to create the observed (3. Moreover, 
t-quarks could decay in the course of the p.t. creating baryon asymmetry on 
the way. 
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Thus we see that if some of t-quark decay channels, or decay channels of 
the produced particles, can accommodate baryon number violation as well as 
CP nonconservation, we have fulfilled all the three basic Sakharov's require- 
ments for a successful baryogenesis. In our model, allowing for B-violation 
coupling in the superpotential, we provide the possibility for the asymmetry 
to be generated, if at least one superquark is lighter than the top quark, so 
that some B-violating channels are kinematically allowed. It is also possible 
to account for a baryon asymmetry in the situation in which the above- 
mentioned superquarks are (slightly) heavier than the top quark, but in this 
case a more detailed analysis of the mechanism through which the superpart- 
ners acquire their masses is needed. Let us remind that baryogenesis in the 
minimal standard model happens to be unsuccessful, in particular, because 
of very weak CP- violation in the MSM, see e.g. the lectures ^3]. I n a SUSY 
extended model CP-violation is not restricted by the CKM matrix and can 
be much stronger. We neglect electroweak nonconservation of baryons in- 
duced by sphalerons because by construction baryogenesis proceeded below 
the EW p.t. 

In what follows we retread in detail every step outlined above. We first 
consider the presence of the baryon- number- violating term in the superpoten- 
tial from a theoretical point of view, and discuss some of its phenomenological 
consequences. Then we will write down the trees of the essential processes, 
and estimate the baryon asymmetry in each case, in a purely particle physics 
way (no plasma, no expansion). Afterward, we will set these estimates in a 
real cosmological environment showing how the necessity of a lower Plank 
mass arises, and find that the asymmetry previously estimated is indeed close 
to the actual one. Finally, keeping in mind restrictions on the model param- 
eters necessary for successful generation of the baryon asymmetry, we will 
discuss the predictions of the model for B-nonconserving processes in particle 
physics, i.e. neutron-antineutron oscillations and heavy quark decays, and 
conclude. 

3 SUSY background 

The 7£-parity violating superpotential, consistent with the gauge symme- 
tries of the SM, as well as with global supersymmetry, made up only of MSSM 
fields, can be written as 

W n = fJL%H u + ^\ ijk LiLjE% + X Hjk LiQjDl + ±\' H *U?D?D° k , (6) 

where H u , L, E c , Q, U c , and D c are the up-type Higgs, left- and right-handed 
lepton, left- and (up- and down-) right-handed quark chiral supermultiplets; 
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the superscript c denotes charge conjugation, and we have suppressed all 
fermion, colour and gauge indices, see e.g. [THITH]. We notice at a glance 
that these couplings violate either baryon or lepton number separately, in 
particular, the first three terms proportional to fi, A and A' have AL = 1, 
while the fourth with coupling A" has AB = — 1, and is the one which we 
will later focus our attention on. For applications of lepton number violating 
terms within the context of leptogenesis see ref. [T^j. The A" term in the 
superpotential (jHJ) leads to the following explicit interaction Lagrangian 2 : 

dnt = -- X" ijk (u* dj d\ + d* k in d) + d* Ui d^j + h.c. (7) 

where the * denotes the complex conjugation of a scalar field, and h.c. means 
hermitian conjugate, and where again fermion, colour and gauge indices are 
suppressed. It is essential for the future estimates that A-' jfc is antisymmetric 
with respect to the last two indices j and k, because of the antisymmetry 
with respect to colour. 

7?.-parity is known to be defined as: 

n = (-if = (-if (B - L)+2s (8) 

where B, L and s are baryon number, lepton number, and spin, respectively. 
It is straightforward to see that if 7^-parity is exactly conserved, then every 
term in W-jz, eq. ©, is strictly forbidden. However, it is not necessary from 
a phenomenological point of view to forbid simultaneously all the couplings 
in ©, since, for instance, the disastrously efficient tree- level diagrams for 
proton decay are not generated if only the A" term is present. Moreover in 
such a theory proton must be absolutely stable. Indeed in this case the lepton 
number is strictly conserved and proton simply does not have any channel 
to decay into 3 . On the other hand, neutron- antineutron oscillations and 
other processes with variation of baryonic number by two, e.g. nuclei decays, 
may become observable with a mild increase of the present-day experimental 
accuracy. If in addition to A" other couplings in eqs. ()6|7j) are also non- 
vanishing they would induce jointly proton decay but one can always take 
them small enough to avoid conflicts with experimental bounds. 

The discussion of all the different combinations of vertices, with their ex- 
perimental consequences as well as the parities associated to them and the 

2 Note that this is not complete since we must include soft SUSY breaking terms con- 
sistent with the above mentioned symmetries, see ref. |15) . 

3 In higher orders of perturbation theory proton may decay due to a possible small 
Majorana mass of neutrino or due to exchange of the heavy Majorana fermion responsible 
for see-saw mechanism, if they exist. 
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underlying theories which provide these symmetries, is beyond the scope of 
this paper, and we address the interested reader to the comprehensive re- 
view [03, in particular to section 2.7, references therein, and ref. \Fj\. The 
main point which we mention here is that we can build a phenomenologically 
acceptable SUSY model with a superpotential which includes the A" term, 
and forbids all the AL ^ ones. These models are usually embedded in (su- 
per)string theories, in which the way the extra dimensions are compactified, 
or, in the case of braneworld, the internal symmetries between the branes, 
provide the necessary symmetries for this configuration to be realised. 

The next point which deserves discussion is the issue of the mass spectrum 
of the MSSM. Since none of the superpartners of the known SM particles are 
observed experimentally, we have only lower limits on their masses which can 
be found in ref. [2] . Theoretically there is no certainty about the form of the 
spectrum, but it is widely believed that one of the two mass-eigenstates of 
s-tops and/or one out of the two of s-bottoms are lighter than all the other 
s-quarks, and that some neutralinos and/or charginos are likely to be even 
lighter. 

Lacking precise knowledge in this sector, we assume as an exemplificative 
toy model a given spectrum compatible with experimental data, and with 
some general requirements on SUSY models. To be more precise, we assume 
either that there is one out of the two mixtures of s-tops, or one of the two 
mixtures of the s-bottoms, which is lighter than the ordinary top quark, 
while all others s-quarks are heavier. We will specify later, when we consider 
quantitatively some concrete model, where the neutralinos and charginos 
reside in this scenario, since different possibilities are allowed and lead to 
(slightly) different results. 

As a final remark let us note that if we take the lower limits, given by the 
Particle Data Group 0, for the superpartner masses 4 , whatever couple of 
(different) particles we choose, their total mass, mi + m?, is either very close 
to or mostly above, the top quark mass. Actually, the only combinations 
with mi + ni2 < m t we can find are Xii or X25 an d one light s-quark or 
chargino 5 . Thus we should not consider decay channels of top quark and 
daughter particles where more than one s-particle are produced, with the 
exceptions mentioned above, since it appears unlikely that these channel are 
kinematically allowed. 

4 Actually these limits have been found assuming 7?.-parity conservation, but we believe 
that possible 7?.-parity non-conserving decays do not change these bounds significantly. 

5 Throughout the paper we denote super-quarks as s-quarks, which should not be con- 
fused with the strange MSM quark. 
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4 Particle physics models and possible decay 
channels 



We outline here four different possibilities of heavy particle decay channels 
determined by a different choice of their mass relations. Below we draw a 
schematic tree for the processes in every situation, and compute the baryon- 
antibaryon asymmetry generated by these decays assuming freely decaying 
particles, thus neglecting the universe expansion and the thermal bath in 
which particles are immersed, where, given the latter, we are neglecting in- 
verse reaction too. At this stage we assume, moreover, that SUSY particles 
are produced only by decays of heavy SM particles and their daughter pro- 
cesses. In the realistic cosmological situation when the i-quark became out 
of equilibrium after the EW p.t. because of m t > T, as is described in 
sec. the inverse decays induced by the superpartners already present in 
the plasma would not be essential because of the Boltzmann barrier and 
thus the assumption of absence of s-particles in the initial state is justified. 

To be more specific, after the EW p.t. there would be heavy t-quarks 
with m t > T but with a non-suppressed number density. There could be 
some other massive SUSY partners as well, also with m s > T. If m s > m t , 
then their decays could contribute to the baryon asymmetry. There are also 
some lighter particles, into which t-quark may decay. The inverse decay and 
resonance reactions with these particles may diminish the asymmetry. How- 
ever, since light particles thermalise quickly, their number density would be 
close to that of massless particles, i.e. n ~ T 3 , and their spectrum would 
be the usual equilibrium one with this T. If this is the case, though it is 
not necessary so, then the inverse decay and resonance scattering of light 
particles would not be important, because the probabilities of these pro- 
cesses are small, ~ exp(— m t /T), while the number density of t-quarks is 
much larger than (m t T) 3 ^ 2 exp(— mt/T) Hence, for every decay process, we 
neglect additional (inverse and resonance scattering) reactions while consid- 
ering the cosmological setup. A more accurate study of the complete kinetics 
of B-nonconservation in the primeval plasma may possibly change the simple 
estimates presented below but not by a large factor. 

In the next section we estimate the baryon asymmetry generated by t- 
quark decays and show that the results are essentially the same as presented 
in subsections 14.11 and 14.21 once we set the Plank mass in a certain low-scale 
range. 
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Figure 1: Schematic representation of the processes involved. As stated in footnoteEl 'u' 
and 'd' stand for generic up- and down-type quarks, while 't' and 'b' are the third family 
top and bottom quarks. Note however that not every combination of quarks is allowed, 
see the discussion in section 0D 

4.1 Lightest s-top 

Let us first consider the scenario in which the only SUSY particles lighter 
than the top quark are the s-top, t±, and the lightest neutralino, x®, where 
the subscript means that we are dealing with only one mixture (labeled 1) of 
a certain kind of s-particles, which is supposed to be the lightest. We begin 
with an equal initial number of top quarks and antiquarks alone, and study 
the evolution of this system. If + m^o < m t , then the possible processes' 
tree is shown in fig. 



10 



We parametrise the branching ratios in the following way 6 : 

r (t -»• W + + d) /Tf = l-x 
r (t -> W~ + d) /Y\ ot = l-x 

r(t-ti + x?)/rf = x 



r (tj -> d + d) / rf = 1 
r f ^ ^d + d) /vf = 1 



r(x?^« + d + d)/rf 



V 



T(x\-u + d + d) /Tf 



2 

1 + Ci 



The last four processes do not conserve baryonic charge and, in the absence 
of 7£-parity breaking interaction (jHJ, t\ and xi would be stable. The dif- 
ference between x and x is induced by CP-violation and is small because 
CP-breaking can manifest itself only in higher orders of perturbation theory 
due to rescattering in the final state, see e.g. [T3ll21| . The same result applies 
to €y., which is non-zero but small. 

Now, if we assume that these particles decay freely and inverse processes 
are negligible, we find that the final asymmetry generated after several char- 
acteristic lifetimes of the decaying particles, is given by summing all the 
quarks and antiquarks in the final states, with their respective baryonic 
charges, multiplied by the branching ratio of the channel into which they 
have been produced. This gives: 

AB = (x — x) — (x + x) ~ e t (9) 

where et = x — x, and in the last equality we have retained only the first 
order terms, since we expect well as to be small. 

If no CP violation is allowed in the standard model particle decay, t — > 
W + q, which means that x = x, we obtain: 

AB = -2xe^, (10) 



6 Hcrc and in what follows 'u' and 'd', together with their antiparticles and super- 
partners, mean up- and down-type quarks of arbitrary family, respectively, while 't' and 
'b' stand for top and bottom quarks of the third family; again the definition extends to 
antiparticles and superpartners. 
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Figure 2: Schematic representation of the processes involved. As stated in footnoteEl 'u' 
and 'd' stand for generic up- and down-type quarks, while 't' and 'b' are the third family 
top and bottom quarks. Note however that not every combination of quarks is allowed, 
see the discussion in sectional The line going from b to b\o is crossed to indicate that 
this channel may be closed if the mass difference does not allow it. 



that is, the asymmetry is entirely generated in the second step 7 , of the CP- 
violating decays of neutralinos. 

Concluding this subsection we notice that if the condition +m.^o < m t 
is no longer satisfied because of a heavier either neutralino or s-top, then the 
baryon asymmetry cannot be generated in this version of the scenario. 

4.2 Lightest s-bottom 

We turn now to the case where the lightest s-quark is the s-bottom b±, and 
the mass spectrum goes as: 

m^o < m~ bi < m t < (11) 

Let us assume again that there are some equal initial numbers of top and 
antitop quarks and no other particles are present. In this case the baryon 
asymmetry can be generated in three steps by the decay of t-quarks and their 
decay products', as sketched in fig. El The last step is actually realised only if 
the chain of inequalities (jllj) is true all the way down. In fact, if the neutralino 

7 It is easy to sec that even if we have more than one channel with AB ^ 0, as in fact 
we have, no asymmetry develops in the primary decay of the top quark, since all these 
B- violating channels have the same AB. 
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is heavier than the b\, the tree of decays changes, because the decay channel 
bi — > b+x® and its charge conjugated one are no longer kinematically allowed, 
and then this process, as well as the subsequent neutralino decay, must be 
cut from the tree, leaving only the AB ^ channel. 

Then, as in the previous case, the branching ratios can be written as 
follows: 



r(t 
r(f 

Tit 



w 
w 



bi + d 



Flt^bt + d 



+ d)/Y 
+ d)/T 

/r 



r (h^b + x 



r(&i 



u 
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/r 
/r 



tot 
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y 
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/r 



tot 

I 



tot 

b 



T(xl^u + d + d)/Y 



tot 

X 



r(xS 



u 



d + d) /r 



/of 

X 



2 

1 + 



Here again the last four processes violate baryonic number conservation but 
in addition the processes 3 and 4 also break it. We see that, with respect 
to the previously considered scenario, here all the three steps can generate 
an asymmetry, since all of them include two different decay channels with 
different AB, either AB c h.\ = and AB c h.2 = ±1 (first and second steps), 
or AB c h.i = 1 and AB c h_2 = — 1 (third step). 

We have chosen this parametrisation since we expect the ratio between 
b\ — > u+d and b\ — > b+Xi t° De small 8 , which in turn means that y and y are 
small, and this gives us a convenient way to express the baryon asymmetry 
generated as: 

AB = (x — x) — (x y — x y) — (x + x — (x y + x y)) — e t , (12) 
where again we could focus on the to obtain: 



AB = -x {{y - y) + (y + y - 2)) ~ -x (e 5 - 2 e 



X) 



(13) 



8 The estimates of the decay rates and the discussion on that will be given in the next 
section. 
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with e- b being y — y, and, as in the previous subsection, the last equalities in 
(jl2j) and (|13[) follow, once one discards higher-order terms, given the small- 
ness of x, x, y, y and e^. 

Actually, as we see in what follows, it is not mandatory that y and y are 
small. Indeed, if this is not the case, we might want to change slightly our 
parametrisation of the branching ratios, and redefine y and y as 1 — z, 1 — z, 
respectively, where now z and z are small, and this gives the following result: 

AB = [xz — x z) — ey.{xz + x z) (14) 

or, with no CP violation in the SM decays, given z — z = ey = —e~ b : 

AB — x (z — z) — x e% (z + z) ~ x ey (15) 

Note that if instead the case of a heavy neutralino is considered, we 
must set y — y — 1, or equivalently z = z = 0. Doing so we find that a 
net asymmetry cannot develop, thus if these scenarios are to work we must 
demand (at least) one light neutralino. As a final comment concluding this 
subsection we point out that there are no kinematically allowed processes 
involving charginos. 

5 Cosmological setting 

In this section we analyse the evolution of an initial equal number density of 
top and antitop quarks in cosmological plasma. We begin by discussing initial 
conditions, then we estimate the decay rates for every process, and study the 
cosmological evolution of the baryon asymmetry. Since the two possibilities 
coming from different mass spectra described above have similar outcome, 
we will focus only on the second scenario, namely m^o < < m t < m t ~ r 

As described in section |2J before the EW phase transition all SM particles 
were massless, except for small thermal masses, since the symmetry was 
unbroken. Following the standard calculations, as given for instance in ref. 
|12j . we find the relation which links the characteristic temperature after the 
transition with the Higgs mass 

T 2 p t ~ 1.5 [m\ + (44GeV) 2 ) (16) 

This relation in fact should be corrected in the context of MSSM, because 
possibly also SUSY particles would take part in the EW transition, but the 
effects are not very strong [TH] . 

Since we require the temperature after the transition to be lower than 
the top quark mass, m t ~ 175GeV, we have the upper bound on the Higgs 
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boson mass m# < 135GeV, while with the current experimental lower bound 
win > 114.4GeV, we find a temperature around 150GeV. Within this range 
of Higgs masses, after the EW phase transition the top quarks and antiquaries 
are the only SM particles which are truly non-relativistic (or better to say, 
semi-relativistic), thus they are interesting candidates for the "parents" of 
out-of-equilibrium decay However we could relax this hypothesis allowing for 
a wider range of Higgs masses, but in that case deviations from equilibrium 
are smaller. Nevertheless a reasonable net asymmetry will develop even in 
this case. 

We proceed through the detailed calculations of the requirements for this 
model to work, completing the discussion given in section |21 In order to know 
if the reactions which we are interested in, are in thermal equilibrium or not, 
we must estimate their decay rates and compare them with the expansion 
rate of the universe, given by the Hubble parameter H = a/a, when T rs m, 
where m is the mass of the decaying particle. If Y/H 3> 1, then equilibrium 
would be established, while if 



the decay would be "frozen". 

The second point we must look at is the stability of top quarks and 
antiquarks, since we don't want these particles to decay before the EW tran- 
sition. This is because otherwise, once the transition takes place, all the 
t and t, now massive, would decay almost instantaneously, and, since they 
are non-relativistic, their number density would be exponentially suppressed. 
Hence, no asymmetry will arise. This requirement can be written as 



that is, the characteristic time of the EW phase transition must be much less 
than the mean life of top and antitop quarks. 

We note here that the above stated requirements of out-of-equilibrium 
decay and long lived particle, in standard cosmology result almost the in 
same numerical constraints, since H — 1/(2 1) or H = 2/(3 1) in radiation or 
dust dominated universes, respectively. Then, when we fix the temperature 
T we have also fixed the Hubble parameter H and the time t at the, almost, 
same numerical value. However this is no longer true if we consider e.g. 
brane worlds or loop quantum gravity cosmologies, and, without going that 
far, even an universe dominated by an inflaton in its oscillating regime. In non 
standard cosmologies H 2 ^ p with p being the energy density of the dominant 
component which fills the universe, generally H 2 ~ L(p) and therefore H ^ 
t~ l , see reviews JHH EDj and references therein. In the case of an oscillating 



ry#< i 



(17) 



Ttoi tp.t. ^ 1 



(18) 
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inflaton dominated universe, we have instead that H ~ T 4 , which again gives 
a different relation between time and temperature, which however depends 
on the specific inflationary model, see j22j. In both these two cases we have 
two different numerical constraints for our requirements (|17jl and ()18j1 . but 
in what follows we assume the normal evolution specified at the beginning 
of this paragraph. 

At this point we have formulated the two basic requirements on the decay 
rates, then we list below the rates we will use afterwards. 



647T 77% 



/ - -\ A" 2 

r(t-5i + d) = j^zimt (20) 
r(S 1 -,5 + xS) - 1^4?^ ( 21 ) 

/- \ A" 2 

r(5 1 -« + d) = ^ (22) 

r x + d + d « ^m| ? 23 

where: y is the EW coupling constant, Zj is a function of the i-th SUSY par- 
ticle mixing parameters, rriq is the mass scale of the heavy superquarks, and 
\"j k is the 7£-parity violating coupling constant relative to the B- violating 
part of the superpotential ©. As before, "u" and "d", mean up- and down- 
type quarks, respectively 9 , while t and b stand for top and bottom quarks; 
the definition extends to antiparticles and their (anti)superpartners. 

Let us now turn to consideration of the demands we want to satisfy. We 
will use the following numerical values. The top quark mass is m t = 175GeV, 
the W mass is m^y = 80.4GeV, the heavy s-quarks mass scale is rrig w 300 
GeV, Z? = Zl = 0.5, though those values are unknown and could well 
be one order of magnitude lower. We set the masses of the unknown light 
s-particles to the values m~ bi m 100 GeV and m^o pa 50 GeV. Furthermore, 
from now on we will drop the apices on A", bearing in mind that it refers to 
the baryon number-violating term in the superpotential. 

Up to this point we have not specified the range which the B-violating 
A couplings could span. We have different processes in which almost every 
different possible A coupling is involved. More specifically, in process (|20|) we 

9 Not every combination is allowed, however. For instance in the process t — * b\ + d, 
only d=down and d=strangc are allowed by the symmetries of the coupling matrices, once 
one takes into account colour. 
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can have A3 31 or A3 3 2, the sub-numbers being the family which the quarks 
belong to, where the first number refers to the up-type quark involved, and 
the other two refer to the two down-type quarks, with Xijk = —\kj- Then, 
in process (j22J) the possibilities are An 3, A123, A213 or A223, while neutralino 
decay, (JZEJ), splits into six different processes determined by A112, A113, A123, 
^2 12, A213 or A 2 2 3- Referring again to the review [T^] (sec. 6), we see that 
some of these couplings, are quite strongly constrained, e.g. An 2 < 10~ 6 , 
A113 < 10~ 3 , while the other could be well above unity. The most stringent 
bounds on A^ are obtained from nucleus stability which is destroyed by 
the baryon-antibaryon transformation. It is discussed in the next section in 
connection with neutron-antineutron oscillations. 

The huge dispersion in the values of the 7£-parity -breaking couplings A^ 
may be related to the similar or larger dispersion in the Yukawa couplings of 
Higgs boson to fermions. 

Once the values of the couplings have been chosen, we are able to de- 
termine the ratios between the different channels of every decay process. In 
particular, we must ensure that the SM top decay is the dominant one with 
respect to the 7£-parity violating one, and we need to know which one of the 
two s-bottom decays is the dominant one in order to choose the most con- 
venient parametrisation of branching ratio as described in section |3J Going 
through the calculations, we obtain for the top quark decay: 



0.12A 2 (24) 
7.5A 2 (25) 



Where both A313 and A323 are allowed in the first ratio, and four different 
couplings A113, A123, A213 and A223 are allowed in the second ratio. Almost 
all of them are allowed to be as large as unity. We see that we indeed have, 
for our values of A couplings, that the B- violating top decay is subdominant, 
while the ratio between the two possible decays of the s-bottom could be 
either bigger or smaller than one for different values of A, which means that 
both possibilities described by the equations (|T2^1 and (JUJ), or ((Hj) and (fTHj) . 
are actually possible. 

We consider now the evolution of these processes in the expanding uni- 
verse after the EW phase transition. We have already pointed out that the 
two basic requirements for the viability of this model are the longevity of the 
top quarks in the early universe, and the out-of-equilibrium feature of their 
decay. Let us fix, for illustrative purposes, the temperature after the EW 



T (t -> 61 + d 
T(t->W + d) 
T (bi — ► u + d 
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transition has occurred to 150 GeV, then, if we want the top quark to decay 
at that temperature in a non-equilibrium regime, we must ensure: 



Tjt^W + d) „ ™™ < i (26) 
H Tp t ~ 2 • 10 5 GeV - 1 ; 



that translates into: 



m Pl < 2 • 10 5 GeV (27) 



where H is given by H 2 = 8tt 3 g X J 2 T 2 /90 mpi, with g* « 110 being the 
number of relativistic degrees of freedom at that time. We have substituted 
r (t — > W + d) for the total the total decay width of t-quark, since this is very 
accurately true, and we have begun from checking of inequality because 
this automatically gives to the top and antitop quarks a long-enough life, 
since in this case the condition t p , t . Y (t — > W + d) < 1/2 holds. Equation (|2"7|) 
gives the basic requirement for this model to work, as it has been discussed 
in section |21 

We now can compute all the other interesting Y/H ratios, and conse- 
quently the baryon asymmetry generated by the processes we have consid- 
ered so far. We know that we can approximate the kinetic equations for 
the processes of interest in different ways, which are selected by the value of 
the ratio Y/H m where H m = H(T = m), m being the mass of the decaying 
particle, see refs. jU|22]. Given all the numerical values specified above, we 
obtain: 

r (t -> W + d) mpi 

H mt ~ 3-10 5 GeV 



r (t 6i + d 



H mt 2.5 • 10 6 GeV 



(29) 



H m , 5 • 10 6 GeV 



(30) 



T ^bi — > u + d 



H mr 10 6 GeV 



(3D 



b 



Y(x°i->u + d + d) 2 mpi 

A Tb^GeV (32) 

All the above expressions show that if condition (|2*7j) is satisfied then 

automatically it would set these ratios (much) below unity. This in turn 
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means that the baryon asymmetry, up to some subsequent dilution processes 
through entropy production, is simply given by: 

P « W~ 2 AB (33) 

where AB is given by the formulas of section |3J Whilst every different specific 
value of the couplings A and the other parameters (mixings and so on) gives 
different AB, we notice that there are several appealing possibilities which 
could easily deal with the experimentally observed value of /3, without requir- 
ing unnatural high values of the various CP-violation parameters involved, 
e ~ 10 -5 , even without allowing CP violation in the SM decay t — > W + b. 
In this connection we notice that despite the fact that up to date no Electric 
Dipole Moments have been detected for neutrons or electrons, thus disfa- 
voring maximal CP violating phases in the MSSM, nevertheless the upper 
bounds [2j could be satisfied also for large values of the CP violating phases, 
provided that some cancelation mechanisms are at work, see ref. j2B] and 
references therein. 

As a final remark in this section, let us note that the decay of the neu- 
tralinos, given the expression would take place at temperatures far 

below BBN. However, the above mentioned decay width is no longer correct 
below the scale at which mpi recovers its classical value, which must be be- 
fore the nucleosynthesis, that is, according to [B], around 10 to 100 MeV. 
Once mpi reaches 10 19 GeV the decay of the neutralino proceeds almost in- 
stantaneously (in comparison with the cosmological expansion rate), thus 
preventing possible troubles with late-time entropy production which could 
destroy the successful BBN predictions. Then, since we don't specify how 
the Planck mass goes up to its standard value or, in other words, which is 
the time scale for this process, we don't know whether neutralino decays are 
in equilibrium or not. However this will not change much the results of the 
scenarios outlined above, since in most of them the asymmetry produced by 
their decays contributes only a small fraction to the total asymmetry, except 
for the possibilities given by eq. (fTOj) and 

6 Baryon nonconservation in particle physics 

As we have mentioned above, one of the main goals of this work to present 
a model of baryogenesis which is capable to explain the observed asymmetry 
(|T|) and simultaneously does not contradict very strict lower bound on the 
proton life-time. As we see below in this section, this model with the chosen 
parameters values leads to noticeable, but still below the existing limits, 
neutron-antineutron oscillations and potentially observable nonconservation 
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of baryonic charge in heavy quark decays. According to [2] the characteristic 
time of oscillations for free neutrons is restricted by 

r nn = l/5m nn > 8.6 ■ 10 7 sec, (34) 

while for bound neutrons the limit is slightly stronger: 

T nn > 1.3 • 10 8 sec (35) 

The latter is obtained from the observed stability of different nuclei at the 
level of 10 32 years. If a bound neutron can turn into antineutron inside a 
nucleus the latter would definitely explode. Non-observation of this process 
allowed to put the above quoted limit. The arguments go as follows. The 
admixture of antineutron to neutron in a nuclei is equal to 

sm9 nfi = 5m/ AE (36) 

where 5m is the amplitude of (n — n) -transition in vacuum and AE ~ 100 
MeV is the energy difference of antineutron and neutron in a nuclei. 
The probability of n annihilation in a nucleus is equal to 

Tnn « (TnnVN sin 2 6 nn (37) 

where a nn is the nn annihilation cross section, V is relative velocity of n and 
n and N ~ m 3 ~ (100 MeV) 3 is the number density of nucleons inside a 
nucleus; a nn V ~ l/^ 2 - Assuming, according to the data, that the life-time 
of nuclei is bounded by r(A — »■ A — 2) = l/T nn > 10 32 years we find 



T n n = 1/Sm > (l/AE)y/a nfi VN ^lO 32 years (38) 

Substituting here numerical values indicated above we obtain the result which 
is incidentally quite close to the limit obtained for free neutrons (|35j). 

In the model under consideration the simplest process which could lead to 
baryon-antibaryon transformation is (|23ft . On a slightly more fundamental 
level it is described by a tree diagram of (wd)-transition into b\. Then the 
latter captures another <i-quark and turns into Xv After that the process 
goes in charge conjugated order: x° transforms into b\ and b\ goes into (ud). 

The amplitude of this process is 

A BS = Xll29 l Z *° Z?h {^)\ (39) 
m~ m^o 

bi a.1 

where C is the charge conjugation matrix and (ipCip) 3 is the product of 3 
spinors of initial state quarks and 3 of the final state antiquarks. The process 
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of neutron-antineutron oscillations would be very efficient if udd quarks (and 
antiquarks) belonged to the first generation. However, this is not so and the 
simplest process includes u and d quarks of the first generation and another d 
quark of the second generation, i.e. the strange quark s. Neutron-antineutron 
oscillations must include transformation with change of strangeness by two 
units, AS = 2, and may proceed with additional second order weak interac- 
tions, see ref. |2~5] . 




Graph 1. Fcynman diagram leading to ns transformation or to AA-oscillations. Here 
every quark is labeled with its proper name. 

We can do the same exercise, which is done above for nn-transformation, 
considering e.g. the process of neutron transformation into S-hyperon. A 
stringent upper bound, Am < 10~ 6 — 1CT 7 , was obtained from nuclei sta- 
bility with respect to this process in ref. (see also the review [TH]). A 
similar process would create AA-oscillations. Using the same arguments as 
presented above we find that the life time of AA-oscillations is bounded by 
approximately r AA > 10 6 sec. A weaker bound in comparison with n — n case 
is related to a larger energy difference of S in comparison with n and to a 
smaller annihilation cross-section of riE. 

Assuming the specified above mechanism for baryonic number nonconser- 
vation we can naively estimate the time of nn-oscillation using the following 
chain of processes: neutron transforms into A by first order weak interaction, 
then A goes into A through the process described above and and A "returns" 
to antineutron again by weak interactions. According to that the time of 
nn-oscillations is equal to 



(m A - m n y 
a 2 



(40) 



where m n and m A are the masses of neutron and A respectively and /z nA is 
the matrix element or weak transition of neutron to A. The latter can be 
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expressed through the amplitude of the decay of A — > pir , using PCAC and 
current algebra [2lj] : 

/i nA « A A ^ p7r -/m^, (41) 

where Aa^ P7I — is the amplitude of the indicated decay of A, which can be 
easily found from the life-time of A (we neglected s-p wave problem because 
it is sufficient for an order of magnitude estimate presented here. Thus we 
find: 

T nn ~ 10 12 r AA (42) 

This is by far beyond any existing and future experimental sensitivity. The 
weakness of the (rm)-transformation is related to the necessity of the second 
order strangeness changing weak transition with AS = 2. 

One can avoid the weak processes with AS = 2 allowing diagrams with 
the exchange of superpartner of VT-boson, wino, W, instead of xi> which 
permits the vertex W —> us. Such diagrams lead to the amplitude (uuddds) 
suppressed only by the Cabibbo-like mixing sin 8 C between wino and strange 
super-quark plus u-quark. Hence the amplitude of (mtdddd)-transition is 
suppressed only by first order weak interaction. But the effect is still weak. 
However, even this is not necessary and the weak processes with AS = 1 can 
be also avoided leading to the allowed time of nn-oscillations just above the 
existing bounds. This makes further improvement of the accuracy in search 
for this process quite promising. 

The mechanism which could lead to noticeable nn-oscillations is the fol- 
lowing. Let us consider the Feynman diagram similar to Graph 1 but with 
substitution of zino, Z, instead of neutralino, xi- Now the superpartner of 
s-quark, s, can transform into non-strange quark d emitting Z. This process 
was impossible with Z-boson due to unitarity of the mass matrix of quarks. 
Because of that the transformation matrix Zqq was always proportional to 
unity matrix both in the mass and flavour eigenstate basis. This is not 
necessarily true in the case of non-diagonal transition of quark and its su- 
perpartner. Though the mass matrices of quarks and s-quarks are unitary as 
well, but they can be different for quarks and their superpartners. This allows 
for the strangeness changing processes in "neutral current" interaction with 
an exchange of zino [2Z|. Correspondingly the amplitude (uudddd) would be 
suppressed with respect to (uuddss) only by a mild factor sin 2 9 C , where 9 C is 
Cabibbo-like mixing angle for s-quark mass matrix. In fact this mixing angle 
may be even not as small as the Cabibbo angle in usual strangeness charging 
decays. This opens window for quite strong nn mixing. 

Let us briefly mention baryonic charger nonconservation in heavy quark 
decays induced by interaction (|f)|7jl . These processes are of course discussed 
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in the literature, see e.g. review [T^]. Our only input here is that the con- 
stants Xijki connecting heavy (t, b) generation with lighter ones, cannot be as 
small as e.g. Am < 10~ 6 . Otherwise, the baryogenesis in this model would 
not be efficient enough. Unfortunately we cannot predict how small X^f. are 
allowed to realize still successful baryogenesis because of uncertainty in the 
value of the effective Planck mass and mass spectrum of SUSY particles. We 
could find the range of couplings, masses, and CP-violating phases which are 
necessary for successful baryogenesis to predict possible range of magnitudes 
of B-nonconservation effects in particle physics but this is a complicated task 
which we will postpone. 

7 Conclusion 

We have discussed here a scenario of baryogenesis based on a SUSY model 
with broken 7£-parity in such a way that only baryonic charge is noticeably 
nonconserved, while leptonic charges are either conserved or their noncon- 
servation is very weak. 

In the model presented here the baryon asymmetry of the universe can 
be produced with a sufficiently large magnitude to agree with observations 
but this demands low scale gravity to ensure the necessary deviations from 
thermal equilibrium. 

On the other hand, proton would not decay despite nonconservation of 
the baryon number because the decay demands nonconservation of leptonic 
charge, L, which may be arbitrary weak. In the spirit of the mentioned above 
hierarchy of baryo-nonconserving couplings grossly increasing from light to 
heavy generations we may assume similar picture for lepton nonconserving 
couplings. If it is true, lepton nonconservation would be weak in electronic 
and muonic sector but could be much stronger in tauonic sector. Lepton 
number nonconserving decays of r would be an interesting feature to search 
for. If leptonic charge is noticeably nonconserved only in r-sector, proton 
decay would be suppressed, even with nonconserved leptonic charge. 

Still there are several B-nonconserving effects (with conserved L) which 
may be observable in direct experiments. Neutron-antineutron oscillations 
may have the characteristic time just above the existing upper bounds or 
particles involving heavy quarks could decay into channels with AB = 2. 
Unfortunately the predictions are not accurate and theoretical uncertainty 
may be as large as several (3-4 or more) orders of magnitude due to poor 
knowledge of the coupling constants, A, necessary for baryogenesis. 

If 7£-parity is indeed broken in this way, then the lightest SUSY particle 
(LSP) should be generally unstable and cannot make the cosmological dark 
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matter (DM). We are not completely sure if, in these frameworks, is impos- 
sible to make sufficiently long-lived LSP with the life-time bigger than the 
universe age but have not yet found rigorous arguments either way. If LSP is 
not DM then other candidates as axion, warm sterile neutrinos, primordial 
black holes, etc could be DM. 

There are several problems yet to study. An interesting one is an inves- 
tigation of proton decay in a model with heavy Majorana neutrino which is 
responsible for light neutrino masses. Another thing which could be done is 
a more accurate study of kinetics of generation of baryon asymmetry with 
different choices of the (s)particle mass spectrum. This however, is rather 
an academic problem because we do not know the appropriate parameters 
(couplings and masses) and accurate calculation with unknown numbers are 
not particularly interesting, at least at this stage. 
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